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Abstract 

The forms of the neutral, non-strange pseudoscalar propagator matrix and 
mixed axial current correlator, (0|T(A^yl^)|0), are discussed at next-to- 
leading (one-loop) order in chiral perturbation theory, and the results com- 
pared to those obtained using QCD sum rules. This comparison provides a 
check of the truncations employed in the sum rule treatment of the current 
correlator. Values for the slope of the correlator with q 2 in the two approaches 
are found to differ by more than an order of magnitude and the source of this 

discrepancy is shown to be the incorrect chiral behavior of the sum rule result. 
12.39.Fe, 24.85.+p, 21.30.+y, 21.45,+v 
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In the traditional treatment of charge symmetry breaking (CSB) in the meson-exchange 
framework, contributions to CSB observables arising from isoscalar-isovector meson mixing 
are obtained under the assumption that the strength of meson mixing is g 2 -independent. A 
number of recent papers — PTOf] , however, have demonstrated that this assumption is suspect. 
Although the quantity in question (the off-diagonal element of the matrix propagator) is 
the element of an off-shell Green function and, as such, not in general invariant under 
allowable field redefinitions (in the sense of Haag's theorem), nonetheless, the existence of q 2 - 
dependence for any particular choice of interpolating field throws the standard treatment of 
CSB into question. Of special interest, among the papers mentioned above, is the treatment 
of p—uj and n—rj mixing using QCD sum rules |7]|| since, without having made any apparent 
explicit choice for the meson interpolating fields, the authors claim to extract the leading 
g 2 -dependence of the off-diagonal propagator element. If true, this result would be extremely 
interesting, suggesting that at least the leading g 2 -dependence, was field reparametrization 
independent, and hence to be incorporated in all treatments of CSB. 

In this paper we critically investigate this claim in the pseudoscalar sector (cf. Ref. |5| ) 
and come to two main conclusions. The first is that it is not possible to extract the q 2 - 
dependence of the off-diagonal propagator matrix element from that of the off-diagonal 
element of the axial vector current correlator matrix. That this is true in general is a 
consequence of the well-known behavior of quantum field theories under allowed field redef- 



initions ]TT| : S-matrix elements are unchanged, but off-shell Green functions are not. The 
off-diagonal element of the propagator is not a physical observable; the axial vector current 
correlator is. They cannot, therefore, in general, be related. The second, and more useful, 
conclusion is that chiral perturbation theory (ChPT) provides formal constraints useful for 
investigating the reliability of the trunctations employed in the sum rule treatment of the 
axial vector current correlator. The reason this is true is that ChPT, being constructed 
solely via symmetry arguments, provides the most general possible representation of the 
physics of the pseudoscalar Goldstone bosons realizing the exact symmetries of QCD and 
breaking the approximate chiral symmetries in exactly the way they are broken in QCD. The 



off-diagonal element of the axial current correlator, therefore, has a low-energy expansion in 
terms of current quark masses, momenta and the low-energy constants of the effective chiral 
Lagrangian, the terms of which can be calculated in a reliable and systematic manner. The 
only approximations enter when one truncates this expression to a given order in the chiral 
expansion. However, even if the expansion converges slowly for a given observable (typi- 
cally, expansion to one-loop order for SU(3)l x SU(3)r is sufficient, but a small number of 
observables are known for which this is not the case - for a general discussion see the recent 
review by Ecker [12|] ) the formal dependence on the current quark masses and momenta 



obtained to a given order is a rigorous consequence of QCD. For the case at hand, namely 
the off-diagonal element of the axial current correlator, we will see below that, while the 
leading chiral behavior of the g 2 -independent part is correctly reproduced by the truncations 
employed in Ref. |J , the leading chiral behavior of the g 2 -dependent piece is not. The sum 
rule result for the slope with respect to q 2 of this correlator, which differs numerically from 
the one-loop ChPT expression by more than an order of magnitude, thus cannot be correct. 

The remainder of the paper is organized as follows. We first revisit the attempt to relate 
the off-diagonal propagator and axial current correlator matrix elements, fixing our notation 
in the process. A general expression for the meson pole contributions to the axial current 
correlator is then given in terms of the meson decay constants and two isospin breaking 
parameters describing the couplings of the axial currents and A % to the physical rj and 
7T°, respectively. The calculation of these parameters is then reviewed, the development 
providing, in addition, an explicit realization of the source of error in the treatment of the 
relation of the propagator and correlator matrix elements in Ref. || . Finally, we compare 
the one-loop ChPT result for the off-diagonal element of the axial correlator with that 
obtained via the sum rule analysis. Using the known dependences of the physical meson 
masses on the current quark masses we show that the truncations employed in the sum rule 
treatment of Ref. H remove the leading chiral behavior of the slope with respect to q 2 , and 
hence are unsuitable for use in treating this feature of the correlator. 

We begin with some notation. Let n, rj represent the physical mixed-isospin ir° and rj 



fields, and 7r 3 , 7r 8 the pure I = 1, flavor octet neutral fields. Then one has, in general, two 
mixing angles, 9 n and 9 V , such that, to 0(9^^), 

7T = 7T 3 + 6^718 7T 3 = 7T - 9 n T] 

7] = -9r,7T3 + 7T 8 7T 8 = 9 V 7T + T) . (1) 

(There is, in general, g 2 -dependent mixing so that Q n ^ 9 V .) To 0(9 ntTt ), i.e., to first order 
in isospin breaking, one defines the isospin-breaking parameter, 9(q 2 ), by 

IW) = ijd±x e ^<0|T(^(*K 8 (0))|0> S (g2 _ ^ _ mg) (2) 

where, from Eqns. (1), ^(g 2 ) = g 2 (6^ — #„.) + (m 2 9 n — rri^9 v ) . The axial current correlator, 
11^, is similarly defined via 

Tlf u = ijd'x e^(0|T(4(*)^(0))|0> = lif{q 2 )q,q v + lif {q 2 )g, v (3) 

where A 3 A* are the 3,8 members of the axial current octet A^ = 57^75 In Ref. || 
the authors evaluate this correlator using QCD sum rules, and attempt to determine 9(q 2 ) 
by considering the pseudoscalar pole contributions to the form factor nf 8 . Invoking the 
PCAC relations (OlA^ |7r 3 (g)) = if n q^ and (0\A 8 |vr 8 (g)) = if^q^ and, evaluating Ylf 8 in the 



pole approximation, they write 



III — fnfr) 



9 n 9^ 



q 2 — m 2 q 2 — m 2 



(q 2 — m 2 )(g 2 — m 2 ) 

which would allow the extraction of 9(q 2 ) if nf 8 were known (e. g. from the sum rules 
treatment). The expression (4), however, is incorrect. The n pole term, corresponding to 
the second term in Eqn. (4), arises from taking 

(0|4|7r(g))= (0|4|(7r 3 + M*)(ff)> = iUq, 

(0\Al\7r(q))= (0\Al\(n 3 + 9„n 8 )(q)) = i/Afc ( 5 ) 
and the first term, corresponding to the rj pole, from 
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(0\A*Mq))= (0|4|(tt 8 - e v 7r 3 )(q)) = if nqit 



(0\Al\ V ( q ))= (0|4|(tt 8 - V3)(g)> = -tfAfc • 



(6) 



The second relations in Eqns. (5), (6) are, however, false. Although (0|y4^|7rs) = (0|v4^|7r 3 ) = 
to leading order in the chiral expansion, beyond leading order, both matrix elements are 
non-zero. This is inevitable at some order, given that Cqcd contains a A J = 1 piece. 
The result that the mixed-isospin current matrix element vanishes at leading order is a 
consequence of chiral symmetry, and does not persist beyond this order. Note, moreover, 
that, at leading order, 7r 3 — n 8 mixing is g 2 -independent. Thus, at this order, 6 n = 9 V and 
the g 2 -dependence of 6(q 2 ) vanishes. To determine the q 2 - dependence, one must go beyond 
leading order, but, beyond leading order, the second of relations (5), (6) are not valid. 
The corrected version of the pole approximation to nf 8 follows from the relations 



(0|4|tt> = iUq» (0|4|tt> = iUeiq^ 
(0|4|??> = if v q„ (0|4|77> = -if n e 2 q. 



(?) 



The isospin-breaking parameters e\ i have been computed to next-to-leading order in the chi- 
ral expansion by Gasser and Leutwyler [TJ| (the calculation is reviewed below). In Eqn. (7), 
f n is the physical 7r° decay constant (which differs from that for the charged pions only at 
O ((m d — m u ) 2 ) |T3| , and f v the physical T] decay constant. ChPT implies f v / f n — 1.3 ]T^] . 
To this (next-to-leading) order, e^, f n and f v are all independent of q 2 . From Eqn. (7) one 
obtains, in the pole approximation, 



38 



eifi 



q 2 — m 2 q 2 — m 2 



Q 2 (e 2 f 2 - ej 2 ) + (mfaf 2 - mle 2 f 2 



2\1 



(8) 



(q 2 — m 2 )(q 2 — ml) 

At leading order ei, e 2 , 9 W and 9 V are all equal to 9 , the leading-order 7r 3 — n s mixing angle 
arising from the leading order off-diagonal term in the meson mass-squared matrix, 

V3(m d - m u ) 



9n 



4(m s 



m) 



(9) 



where m = (m u + md)/2. Also at leading order f n = f v = F so that, at this order, the 
numerator in Eqn. (8) reduces to F 2 8q, which is identical to the leading order expression 



for f-Kf n 9{<f)- The expression (4) is thus correct at leading order in the chiral expansion. 
Beyond leading order, however, e 2 f 2 7^ fnf^ri and €\f 2 7^ fnfrj^n (see below), and knowing 
the numerator in Eqn. (8) no longer allows one to obtain that of Eqn. (3). In particular, 
the crucial g 2 -dependence of 0(q 2 ), which arises only at next-to-leading order and beyond, 
is completely inaccessible from the g 2 -dependence of the numerator of (8). 

We now provide explicit expressions for all quantities appearing in the above discussion, 
at one-loop order in the chiral expansion, obtained using the low-energy effective chiral 
Lagrangian, C e ff, of Ref. |13| . All notation is as defined there. We also keep terms only to 
first order in isospin breaking, i.e., drop terms of O ((m<f — m u ) 2 ). 

Let us first discuss the meson propagator matrix, for the standard choice of meson fields 
13] . 713 , 7Tg^ are the bare 7r 3 , ir s fields. The bare propagator matrix, to one-loop order, 



is of the form 



A(g) 



Zz\q 2 -mi) £ 38 



v 1 



38 



(10) 



where Z 3 , Z s and m 2 , m\ are the ir^ , ir^ one-loop wavefunction renormalization constants 
and squared- masses in the isospin symmetry limit. Expressions for m| and wig in terms of the 



parameters of £ e // are given in Ref. [|13|] , and those for Z 3 , Z 8 and £38 are easily computable 
from Ceff. In what follows we write Z 3 X = I+Z3, Z^ 1 = l + z 8 and £ 38 = £(p) +^ (1)0 + £ {1)1 q 2 . 
£(0), and i^ 1 are g 2 -independent and 0(m^ — m u ). £t \ is of leading chiral order and 
z 3) z 8) £^°, i^ 1 next-to-leading order, m 2 and m§ contain both leading and next-to-leading 
pieces. The renormalized 7r 3 , 7r 8 fields, 7r£ and 7Tg, are obtained, as usual, by 



7r; 



7T„ 



(z- i/2 : 



-(0) 
3 

.(0) 



33 "3 + (Z~ 

-w I 1 7-1/2 

83^3 + l Z 



l/2\ _(0) 
^38^8 



, 7T (0) 
'88^8 



(11) 



where Z is the wavefunction renormalization matrix. From Eqn. (10) one has 



/ 



Z -l/2 



1 + ^3 W ] 



1 \ 



v ¥ {1)1 1 



n Z 8 



(12) 



and from this one obtains the renormalized propagator matrix, to 0(md — m u ), 

-i 



A (r) (?) 



I q 2 



VTin 



-m 



-m 



te 



te 



(13) 



where m 2 e 



-[£(0) + + ^^(mj + wig) — \1(q){zz + z&)]. Using the explicit expressions 
for Z3, Zs, £(o), E^ 1 ' and E^ 1 ' 1 , one may verify that m 2 e is now finite, and diagonalize the 
mass-squared matrix to find the physical, renormalized 7r°, 77 fields. Writing 



7T° = < + P < 

rj = -e r Ti r t + K 

with 6 r = —m 2 e /(ml — nij), one finds, explicitly, 



(14) 



9 r — 80 



1 — 3fi n + 2fi K + ji rj + 
32Bo(m s — m) 



'3mg + m 2 
64vr 2 F 2 



1 + 



m; 



log(m^/m^) 



F 2 



-(3L r T + Ll 



= 9 (1 + 56 r 



(15) 



where fip = [m 2 p\og(m 2 p / fj?)] /32tt 2 F 2 with \i the (dimensional regularization) renormaliza- 
tion scale, L r k are the renormalized low-energy constants as defined in Ref. |13| , m 2 K is the 
average K mass-squared and F, B Q are the two constants appearing in the lowest order part 
of the chiral Lagrangian (F = f n in leading order). The expression (14) is not yet of use to 
us since the 7rJ* and 7Tg fields have mixed isospin. We may re- write (14) as 



7T 



(i+2^x4 0) +^n 
:i + ^ 8 )(-W 0) +4 0) ) 



.(Oh 



(16) 



where 
81 = 6> 

+ 

9 2 = 9 



, 7 4 

1 - + -UK + Hr, + 



3m 2 + m 2 



64vr 2 F 2 



1 + 



^,2 



m 



A" 



log(m 2 r /m 2 ^) 



'm 2 — m 2 
v 64vr 2 F 2 
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1 + \og(m 2 K /^ 2 



32(171% — m 



= 00(1 + 66!) 



3 "/V + + Hr, + 



'3m 2 + m 2 



' 2 2 
777^ — m 



64tt 2 F 2 



- (l + lo, 



64tt 2 F 2 / 

m K /n 



1 + 



m; 



m 2 ^ — m 2 



log(mJm K ) 



m 



F 2 



^-{3U 7 + L 



9 {l + 66 2 ). (17) 



The £6*1,2 are next-to-leading order quantities. If £38 had had no g 2 -dependence, we would 
have had 9\ = 62 at this point. To compare with Eqn. (1) we recast the above results in 
terms of the isospin-pure renormalized fields, 713°^ and vrg ^, i.e., the renormalized fields in 
the isospin symmetry limit, tv^ 7 " = (1 + §23)713 and 7Tg°' )r = (1 + §z 8 )7rg°' ) , obtaining 

rt = -9 (1 + 8B 2 + \z, - \z^f r + . (18) 

Thus, formally, we find 6 V = (1 + 6§i + \z z - ~z 8 ), 6 V = 0„(1 + 56 2 + \z s - \z % ). 

Let us turn to the axial current correlator, If^,. Although the low-energy representation 
of the product A^Af, contains contact terms, these do not affect ITf 8 to next-to-leading order, 
so we may restrict our attention to the product of representations of the individual currents. 
These are easily read off from C e ff of Ref. [TB|] , with the result 



+ 7^(47r + 7r- + K+K~ + K^K^d^ 
or 

K = -Fd^f ] -^(m s + 2m) L,d^ - ^(2m s + m) L 5 d^ 



8B (m d - m u ) r Q (0) , 1 



L B d^> + -j=-^(K + K~ - K*K°)d^> + ■■■ (19) 

where the + • • • indicates terms higher order in the fields, and other terms which do not 
contribute to ITf 8 to one-loop order. As claimed earlier, one sees, beyond the leading order 
terms —Fd^n^ in A 3 ^ and —Fd^n^ in A s ^, terms which will couple A 3 ^ to 7Tg and A^ to 7^. 
From the expressions (19) for A 3 ^ 8 and the inverted form of the relations (16), i.e., 

rrf = (1 - i*s)7r - O (1 + 89, - ±z 8 ) V 



7T 



(0) 



1 x 1 



l--z 8 )r / + ^ (l + ^ 2 --^)7r° (20) 
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(valid to 0(md — m u ) , and to next-to-leading order in the chiral expansion) one may easily 

verify the results of Gasser and Leutwyler for f n , f v , e\ and e 2 . The following equivalent 

forms display explicitly the problems with the second of Eqns. (5), (6): 

/n |,3|/u_ ■ r - f (a , i 8B o( m d ~ m u ) {L K + - L K o) \ 
{Q\WVW) = ~ % fv e *% = -i U% K + C12 + L 5 + -^- 2 I (21) 



a8\ Of \\ - ■ f f (n 8B (m d -m u ) (L K + - L K o)\ 

(?)> = « /ttCi^ = I f v % [On ~ C 12 L 5 -^=— I 



(22) 

where Lp is the tadpole loop integral Lp = fx 4 ~ d J -^dj^rz^r^ and C12 = 0q{59\ — §§2 + Z3 — 
z$). The first terms on the right hand sides of Eqns. (21), (22) are the only ones retained 
in the treatment of Ref. || . The remaining terms are present due to the direct couplings 
of ir^ r (the leading piece of 77) to A? , at 0(m d — m u ), and of 713 r (the leading piece of 7r°) 
to A®, at 0(rrid — m u ). As already expected on general principles, the second of Eqns. (5), 
(6) are seen explicitly to be incorrect at next-to-leading order as, in consequence, is Eqn. 
(4), which is based on them: the isospin-breaking parameter, 6(q 2 ), of the pseudoscalar 
propagator matrix cannot be extracted directly from the mixed axial current correlator. 

It is important to stress that this conclusion is inescapable, independent of whether or 
not the chiral expansions for ei^ are well-converged at one-loop order. The g 2 -dependences 
of nf 8 and 6(q 2 ) arise only at one-loop, and not at leading, order in the chiral expansion. 
Eqns. (21), (22) demonstrate that terms of the same order as those kept must be dropped 
for the g 2 -dependent piece of the numerator of Eqn. (8) to reduce to f-Kf^O^q 2 ). Since such a 
procedure would be inconsistent (moreover the additional terms are not small, numerically), 
we confirm that the general argument, based on the behavior of quantum field theories under 
field redefinitions, is not accidentally evaded in the pseudoscalar sector. 

Let us now turn to the more interesting point of using the chiral expansion as a constraint 
on the truncations employed in the sum rule treatment of the axial current correlator. For 
this we will give the explicit expressions for nf s 2 to one-loop in ChPT and compare the 
formal quark mass dependences of these expressions with those of the results of the sum 
rule analysis. The latter may be obtained from the results of Ref. || (written in terms of 
the physical 7r°, 77 masses) by using the known leading- plus next-to- leading expansions for 



m WiV in terms of the quark masses |13] . It is important to remember that the axial current 
correlator is a physical object. As such, its dependence on q 2 is not to be identified with the 
g 2 -dependence of the angle which diagonalizes the inverse propagator matrix (as obtained in 
Refs. ) • That angle, like 9(q 2 ) (though distinct from it), is dependent on the choice 

of meson fields. 

The low-energy representation of the product of axial currents, A^Af,, consists of two 
parts, (1) the product of the low-energy representations of the individual currents, as given 
in Eqns. (19), and (2) a set of contact terms associated with those terms in £ e // quadratic 
in the external axial fields. To next-to-leading order, explicit calculation shows that the 
product of the representations of the individual currents contributes only to nf 8 , and the 
contact terms only to IT^ 8 . For H| 8 one finds 



n 3 8 _ B (m d - m u ) 



l 2 



3 (log(m^//i 2 ) + l) - 8U 5 



32tt 



(23) 



Although not discussed in detail in Ref. || , the pole approximation employed there leads to 
a n^ 8 (as defined here) which is proportional to q 2 and hence vanishes as q 2 — » 0. It follows, 
from Eqn. (23), that Hi 8 cannot be properly modelled using only pseudovector meson pole 
terms, as in Ref. || . Turning to the remaining form factor, Ilf 8 , we note that, to this order 
in the chiral expansion, the pole approximation is exact. The correct one-loop expression 
for II'P is, therefore, that already written down in Eqn. (8), with /„., f v , e% and e 2 the one- 
loop expressions of Ref. [|13j . In obtaining numerical results below we have rescaled the 



values of quoted in Ref. by a factor of 1.22, in order to account for the larger-than- 



expected violations of Dashen's theorem |14|-|1(J . We thus have as input e\ = 1.67 x 10 2 , 



e 2 = 1.35 x 10 2 , f w = 93 MeV, and = 1.3. Then, setting g v = f 2 e 2 and g n = f 2 e x 

and rewriting Eqn. (8) in the form 

n 38 = ( _Jh 9-k \ ( q 2 (9y - 9n) + {9*™% - g«ml) \ 

1 \q 2 — m 2 q 2 — m 2 ) \ (q 2 — m 2 )(q 2 — m 2 ) J 

to facilitate comparison with Ref. , we find g v = 197 MeV 2 , g n = 144 MeV 2 , to be 
compared with the values g v = 143 MeV 2 , g n = 139 MeV 2 extracted via the sum rule 
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analysis || . (These values have been rescaled from those quoted in Ref. || by the same 
factor of 1.22 discussed above.) To leading order in the Borel mass, M B of Ref. || , 

m u (0\uu\0) - m d (0\dd\Q) 
9^ = 9v = 7^r~ 2 ^ (25) 

which, since (0|mt|0) = (0|g?g?|0) = —F 2 B to leading order in the chiral expansion, implies 
(again to leading order) g w = g v = 9 F 2 , the leading order chiral result. The procedure 
of truncating the operator product expansion (OPE) for A 3 ^Af, at operators of dimension 
6 and dropping terms of 0{m 2 ) employed in Ref. || thus correctly reproduces the leading 
order chiral constraints. The interesting g 2 -dependence of U^ 8 , however, enters first only 
at next-to-leading order in the chiral expansion and here, comparing the results of the two 
approaches for the coefficient, g v — g n , of q 2 in the numerator of the last expression for nf 8 
in Eqn. (24), ChPT to one loop predicts g v — g n = 53 MeV 2 , whereas the results of Ref. || , 
rescaled for the violations of Dashen's theorem, produce g rj — g % = 3.9 MeV 2 , more than an 
order of magnitude smaller. If we consider the leading chiral behavior implicit in the sum 
rule result we can easily discover the source of the discrepancy. Using the expressions from 
Ref. H to evaluate g v — g^, one obtains, to 0{md — m u ), 

* - fc = <0|,- 5 |0> (M±^> + 0(1/K)) (26) 

where (0|gg|0) = (0|uw|0) ~ (0\dd\0), Mb is the Borel mass and the higher order terms 
are also higher order in the meson masses (and hence higher order in the chiral expansion). 
Using the known chiral expansions of the meson masses, one obtains, for the leading term 
in the chiral expansion of g v — g n , 

n n -a F 2(% B 2 {m s - m)(m s + 2m) \ 

where + ■ • • represents terms higher order in the chiral expansion. In contrast, making the 
same expansion for g v = f 2 e 2 and g n = j 2 e\ using the relevant expressions from Ref. JT3] , 
one finds 

9r,-9«- VqF y 87i2F2 iog(m K /fi ) — + — L 5 H J ■ (28) 
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The result of Eqn. (27) begins at 0(mV), that of Eqn. (28) at 0{m q ) (plus the ubiquitous 
chiral logs): the truncation scheme of Ref. [Q has thus removed the leading contribution to 
the slope of Hf 8 , which accounts for the smallness of the numerical result. It would be of 
interest to determine which aspect of the truncation scheme is the source of the problem. 
The situation here appears analogous to that encountered in the sum rules analysis of the 
nucleon mass, where the usual approximate treatment of the continuum contributions to 
the phenomenological side of the sum rule leads to violations of the known chiral behavior 
of Mn with m, and only a more careful analysis restores the correct next-to-leading order 



behavior 117,18 



In summary, we have demonstrated (1) that the mixing parameter, 6(q 2 ), of the pseu- 
doscalar meson propagator matrix cannot be extracted, except at leading order in the chiral 
expansion (where it is independent of q 2 ), from a treatment of the mixed axial correlator 
11^, and (2) that the QCD sum rules treatment of 11]^, in the approximation of dropping 
terms of 0(m 2 ) and truncating the OPE to operators of dimension 6 or less, fails completely 
at next-to-leading order. Since the pole approximation for nf 8 breaks down beyond next-to- 
leading order, the sum rule method appears unlikely to provide any information about the 
g 2 -dependence of Ifi® beyond that already known, at next-to-leading order, from the chiral 
expansion, unless the multi-meson continuum is taken into account in the phenomenological 
side of the sum rule analysis. One could, however, profitably use the results of ChPT to 
provide useful constraints on the application of the sum rules method to few-body systems. 
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